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■ Abstract. Coronal structure and coronal heating are intimately related 
CS| ! in magnetically active stars. Coronal structure is commonly inferred from 

Q ■ radio interferometry and from eclipse and rotational modulation studies. 

O , We discuss to what extent flares may be responsible for coronal structure 

Q ' and global observable properties in magnetically active stars. 

in ■ 

T— I \ 

^ \ 1. Introduction 

OO ' 

\ What are stellar coronae? The answer to this simple-minded question is deeply 

. rooted in the interplay between magnetic field generation and energy release 

CN \ in stellar atmospheres. The nature of this interplay is surprisingly difficult to 

■ grasp. While observables such as coronal plasma temperatures, emission mea- 
\ sures, electron densities, non-thermal particle densities and in cases magnetic 

field strengths have painstakingly been measured, monitored, and correlated, 
our question remains one of understanding cause and effect in a complicated 
array of interactions. There is little doubt that magnetic fields not only confine 
O \ coronal plasma but also heat this plasma once they interact due to shearing and 

' twisting as a consequence of surface convective motions. Coronal heating thus 

\ intimately depends on magnetic structure; conversely, magnetic reconnection 

• determines the geometry of radiating magnetic structures. 

' 2. Coronal Structure: Magnetospheres, Flares, Active Regions 



2.1. Radio Interferometric Imaging of Stellar Coronae 

Radio Very Long Baseline Interferometry (VLBI) has produced some of the best 
evidence of stellar coronal structure related to non-thermal electrons around 
active stars, on milliarcsecond (mas) scales. Since the electrons are trapped in 
magnetic fields, they outline their overall structure. The achievable resolution is 
still only of order of a stellar diameter for nearby active stars or of the intrabinary 
distance for the closest RS CVn and Algol binary systems. 

RS CVn and Algol-type binaries often reveal a bi-modal compact core + 
extended halo radio structure, the halo being of a total size comparable to 
the binary system size (e.g., Mutel et al. 1985). The extended structure is 
resolved in two oppositely polarized radio lobes in Algol that are aligned with 
the rotation axis (Mutel et al. 1998; Fig. lb). The large-scale structure seems 
to be based on globally ordered magnetic fields as suggested by polarization 
gradients (Beasley & Giidel 2000). Other stellar classes exhibiting large-scale 
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Figure 1. Radio coronal structure: Top left: VLBI image of UV Get 
(resolution ~ 0.8 mas; after Benz et al. 1998). Top right: VLBI image 
of Algol (Mutel et al. 1998; units: mas). Bottom left: Model for a 
radio magnetosphere of a young star (Andre et al. 1988). Bottom 
right: Magnetospheric model after Mutel et al. (1998). 

magnetospheres include pre-main sequence and magnetic chemically peculiar 
Bp/Ap stars (Phillips & Lestrade 1988; Philhps, Lonsdale, &: Feigelson 1991, 
1993; Andre et al. 1991, 1992; Smith et al. 2003) with sizes up to lOR^. 

The spatial and polarimetric boundary conditions are satisfied by global 
dipole-like magnetospheres somewhat resembling the Earth's Van Allen belts 
(e.g., Andre et al. 1988; Morris, Mutel, & Su 1990; Linsky, Drake, & Bastian 
1992; Figure Ic). Stellar winds escaping at higher latitudes could draw the 
magnetic field lines into a current sheet configuration in the equatorial plane. 
Particles accelerated in that region travel back and are trapped in the equatorial 
magnetospheric cavity. Additional, dynamic ingredients may be required: The 
compact cores sometimes seen during enhanced flux episodes are likely to be ex- 
panding flares that evolve into extended halos where the low-density population 
of accelerated electrons slowly cool on time scales of days (Mutel et al. 1985). 

VLBI techniques have been more demanding for dwarf stars, but spatially 
resolved observations of some M dwarfs also reveal coronae of a size much larger 
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than the stellar diameter (Alef, Benz, Sz Giidcl 1997; Pestalozzi ct al. 2000). For 
example, the dMe star UV Get B is surrounded by a pair of giant synchrotron 
lobes, with sizes up to 2.4 x 10-^° cm and a separation of 4—5 stellar radii along the 
putative rotation axis of the star, suggesting very extended magnetic structures 
above the magnetic poles (Benz, Conway, & Giidel 1998; Fig. la). In this case, 
trapped high-energy electrons radiate most efficiently in the converging magnetic 
fields over the polar regions (see also Mutel et al. 1998 for a similar model for 
the magnetic field of Algol, Fig. Id). 

The prime contribution of VLBI to coronal physics has thus been in the 
recognition of extended, globally structured magnetic fields containing a variable 
number of high-energy electrons. The origin of the latter is not definitively 
known, and our knowledge of the connectivity of the magnetospheres to the 
surface magnetic fields is model-dependent. 

2.2. Eclipses 

The most productive methods of coronal structure mapping are based on (usu- 
ally non-unique) image reconstructions from eclipse light curves (White et al. 
1990 = W90; Siarkowsi 1992; Schmitt 1996; Giidel et al. 2003a). 

Studies of extremely active binaries support a bi-modal thermal and ge- 
ometric structure, somewhat reminiscent of structure inferred from VLBI and 
clearly unlike the non-flaring solar corona (Walter et al. 1983 = W83; White et 
al. 1986 = W86; Stern ct al. 1992; Ottmann 1993, 1994; Ottmann, Schmitt, & 
Kiirster 1993 = 093): i) Judged from rather shallow or absent X-ray eclipses in 
some RS GVn or Algol-type binaries, a significant fraction of the magnetically 
enclosed X-ray emitting plasma must be extended on scale heights ^ i?*; these 
structures are often associated with subgiant components (e.g., in AR Lac or 
Algol). They seem to be predominantly very hot in X-rays (» 10 MK), remi- 
niscent of flaring plasma. The non-detection of radio eclipses (Doiron & Mutel 
1984 for AR Lac; van den Cord et al. 1989 and Mutel et al. 1998 for Algol; 
Alef et al. 1997 for YY Gem) similarly suggests an extended magnetosphere 
also detected and resolved by VLBI (§2.1). ii) A second type of X-ray structure 
found in both evolved and main-sequence binaries shows high pressure (up to 
several 100 dyne cm~^) and small, "solax-like" scale heights. Although they may 
be cooler than the extended structure, their pressure is reminiscent of flares. 

The emission is inhomogeneous, concentrated in active regions in particu- 
lar on the leading hemispheres of the stars (W90, 093) or between the stars 
(Culhane et al. 1990; Siarkowski 1992; Siarkowski et al. 1996; Pres, Siarkowski, 
& Sylwester 1995) and at low to mid latitudes (W83, W86; W90; 093; Giidel 
et al. 2001a = GOla; Fig. 2e), suggesting that the "active" filling factors are 
modest (up to 10 — 20%). The possibility that some of the material links the 
two stars (Fig. 2c) opens up speculations on magnetic reconnection between the 
two coronae (see, e.g., Uchida & Sakurai 1983). Such has also been suggested 
from radio eclipse observations of the Algol-type binary V505 Sgr (Gunn et al. 
1999) and the pre-cataclysmic binary V471 Tau (Lim, White, & Cully 1996). 

It is not a priori clear in which type of structure large flares occur. For 
M dwarfs, hydrodynamic simulations of large flares suggest active region sizes 
that cover a substantial fraction of the stellar surface (Reale ct al. 2003). VLBI 
results suggest that flaring evolves from compact to extended structures, and this 
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Figure 2. Stellar X-ray eclipses and map reconstructions: (a) YY 
Gem eclipse (GOla); (b) a CrB eclipse (Giidel et al. 2003a). Bottom, 
clockwise from upper left: AR Lac corona (Siarkowski et al. 1996); a 
CrB corona (Giidel et al. 2003a); solar image for comparison (Yohkoh 
mission of ISAS, Japan); YY Gem corona (phase 0.875, GOla). 
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seems to be supported by eclipse observations of flares on VW Cep and Algol 
that indicate relatively compact emission sites with sizes of « 0.85R^ (near one 
pole, Choi &: Dotani 1998), ^ 0.6i?* (scale height, near one pole, Schmitt Sz 
Favata 1999) and O.IR^ (Schmitt, Ness, & Franco 2003). Densities of 5 x 10^° - 
several times 10^"^ cm~^ are found, similar to solar flares. 

While the above examples are all extremely active stars, a transition to 
more solar-like behavior is expected as the activity declines. A particularly fa- 
vorable configuration is ofi"ered by a CrB B, an intermediately old (~ 500 Myr) 
solar analog that is totally eclipsed by an X-ray dark A star. Image reconstruc- 
tions from light curves indicate a patchy structure with compact active regions, 
revealing electron densities up to a few times 10^*^ cm~^ but no evidence for an 
extended corona (Schmitt & Kiirster 1993; Giidel et al. 2003a; Fig. 2). 

2.3. Rotational Modulation 

In the absence of eclipses, a more general approach to coronal structure is pro- 
vided by monitoring of rotational modulation. The Sun often shows rather 
pronounced X-ray rotational modulation as a small number of active regions 
alternately rotate into and out of view. However, rotational modulation is ex- 
ceptional among stars, one of the main reasons being that most of the brightest 
examples are highly active; such stars are probably covered with many active 
regions, and intense flaring often covers up low-amplitude modulations. Such 
problems have hampered identification of periodic X-ray modulation in the pro- 
totypical young active star, AB Dor (White et al. 1996 and references therein; 
Giidel et al. 2001b; but see Kiirster et al. 1997), although radio observations 
reveal two emission peaks per rotation seen repeatedly over long time intervals. 
They probably relate to preferred active longitudes (Lim et al. 1992). 

Somewhat less active stars may be better-suited targets. EK Dra, a young, 
active solar analog, shows rotational modulation both in X-rays and radio (Giidel 
et al. 1995). The depth and length of the modulation constrains the X-ray 
coronal height and the electron densities (> 3 x lO^'' cm~^) and demonstrates 
that much of the emitting material is concentrated in extended "active regions" . 

The X-ray luminosity is also below its possible maximum for stars in the 
"supersaturation regime" where the X-ray luminosity decreases with increasing 
rotation velocity. The cause is not clear, but a deep modulation of the super- 
saturated young star VXR45 suggests that extreme activity in these stars is not 
due to complete coverage of the surface with active regions (Marino et al. 2003). 

Among evolved stars, the RS CVn binary HR 1099 (Kl IV -|- G5 V) has con- 
sistently shown X-ray/EUV rotational modulation, with a maximum at phases 
when the larger K subgiant is in front (Agrawal & Vaidya 1988; Drake et al. 
1994; Audard, Giidel, & Mewe 2001). The rotationally modulated material can 
therefore be located - in contrast to results from eclipsing studies described above 
- on the K star hemisphere that faces away from the companion (Audard et al. 
2001, using information from Ayres et al. 2001, see §2.4). 

2.4. Doppler Information 

Doppler information from X-ray spectral lines may open up new ways of imaging 
stellar coronae as they rotate, or as they orbit around the center of gravity in 
binaries. First attempts are encouraging although the instrumental limitations 
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are severe. Ayres et al. (2001) find Dopplcr shifts in lines of HR 1099 that 
clearly agree with the line-of-sight velocity of the subgiant K star. Periodic 
line broadening in YY Gem, on the other hand, suggests that both components 
are similarly X-ray luminous (GOla). Huenemoerder et al. (2003) find Doppler 
motions in AR Lac to be compatible with coronae on both companions, close to 
the photospheric level. For the contact binary 44i Boo, Brickhouse et al. (2001) 
find periodic line shifts, and from the amplitudes and rotational modulation 
(Brickhouse &: Duprcc 1998), they conclude that two dominant X-ray sources 
are present, one very compact and one extended, both located close to the stellar 
pole on the larger companion (see also Brickhouse, these proceedings). 

3. Energy Release and Coronal Structure 

The study of coronal structure confronts us with several problems that are diffi- 
cult to explain by scaling (up or down) of solar coronal structure: i) Character- 
istic coronal temperatures increase with increasing magnetic activity (Schrijver, 
Mewe, & Walter 1984; Schmitt et al. 1990; Giidel et al. 1997). ii) Characteristic 
coronal densities increase analogously, at least in main-sequence stars (Giidel et 
al. 2001b; Ness et al. 2001, 2002). iii) The maximum X-ray luminosity exceeds 
the level expected from complete coverage of the surface with solar-like active 
regions by up to an order of magnitude (Vaiana & Rosncr 1978). iv) Radio ob- 
servations reveal persistent high-energy non-thermal electrons in magnetically 
active stars even if the lifetime of such a population should only be tens of min- 
utes to about an hour under coronal conditions (Giidel 2002). Several of these 
features are reminiscent of flaring, as are some structural elements in stellar 
coronae (§2). The remainder of this paper is devoted to the hypothesis that 
coronal structure is intimately related to explosive energy release, i.e., flares. 

3.1. Correlations between Coronal Emissions 

The time-averaged power from optical U-band flares correlates linearly with 
the low-level, "quiescent" X-ray luminosity (Doyle & Butler 1985; Skumanich 
1985). The latter itself correlates with the rate of X-ray flares exceeding a given 
energy threshold (Audard et al. 2000). Such would be expected if quiescent 
emission were physically related to the overall flaring rate in the active corona - 
for example if the detected X-ray flares represented the "tip of the iceberg" , the 
iceberg being a distribution of flares of various amplitudes, the ocean being the 
combined emission of a large number of superimposed small flares. 

Giidel & Benz (1993) discuss a global relation between non-thermal radio 
luminosities of active stars and their "quiescent" X-ray luminosities. Since the 
lifetime of MeV electrons in coronal magnetic fields requires frequent accelera- 
tion, a possible explanation again involves stochastic flares: the flare-accelerated 
electrons could themselves act as the heating agents via chromospheric evapo- 
ration. The smoking gun comes with the observation that the total radio and 
X-ray outputs of solar flares follow the same correlation (Benz & Giidel 1994). 

Since i) flares heat plasma, ii) increase the electron densities, therefore iii) 
increase the luminosity per surface area, and iv) amply produce non-thermal 
electrons, could it be that the observed phenomenology and structure on mag- 
netically active stars is due to an ensemble of stochastic flares? 
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3.2. Coronal Variability: Quiescent versus Flares? 

Suggestive evidence has been gathered from hght curves. A strong correlation 
between H7 flare flux and simultaneous low-level X-ray flux in dMe stars sug- 
gests that a large number of flare-like events may be responsible for the low-level 
emission (Butler et al. 1986). At least in main-scqucncc stars, there is little evi- 
dence for constant emission in high-sensitivity X-ray data (Audard et al. 2003). 
Frequent faint, flare- like X-ray fluctuations (Fig. 3b) are often accompanied - in 
fact slightly preceded - by optical (U band) bursts, the latter being a signature of 
the initial bombardment of the chromosphere by high-energy electrons (Giidel et 
al. 2002). Stochastic chromospheric evaporation could thus be the process that 
enriches the corona with plasma. In any case, any successful model of coronal 
structure and heating must account for the observed continuous variability. 

3.3. Coronal Flare Statistics 

The suggestion that stochastically occurring flares may be largely responsible 
for coronal heating is known as the "microflare" or "nanoflare" hypothesis in 
solar physics (Parker 1988). Observationally, it is supported by evidence for the 
presence of numerous small-scale flare events occurring in the solar corona at 
any time (e.g., Lin et al. 1984). Their distribution in energy is a power law, 

where dN is the number of flares (per unit time) with a total energy in the inter- 
val [E, E + dE]. If a > 2, then the energy integration (for a given time interval) 
diverges for i^min 0, i.e., by extrapolating the power law to sufficiently small 
flare energies, any energy release power can be attained. This is not the case 
for a < 2. Evidently, then, one needs to measure the energy distribution of a 
statistically relevant number of flares. Solar studies have repeatedly resulted in 
a values of the order of 1.6-1.8 for ordinary solar flares (Crosby, Aschwanden, 
8z Dennis 1993), but recent studies of low- level flaring suggests a = 2.0 — 2.6 
(Krucker & Benz 1998; Parnell & Jupp 2000). 

Relevant stellar studies have been rare. Early results converged to a < 2 for 
stellar samples (Table 1). Mixing stars at different distances and with different 
luminosities can introduce strong statistical bias. To avoid this, Audard, Giidel, 
& Guinan (1999) and Audard et al. (2000) applied a flare search algorithm to 
EUVE light curves of individual active main-sequence stars, taking into account 
flare superpositions and various binning to recognize weak flares. These results 
indicate a predominance of relatively steep power laws including a > 2. 

The identification of weak flares close to the low-level emission level is an 
ill-defined problem. Limited signal-to-noise ratios add to the problem. Some of 
these complications were overcome by fully forward modeling the superposition 
of a statistical ensemble of X-ray or EUV fiares (Kashyap et al. 2002; Giidel 
et al. 2003b) including an analytical formulation (Arzner Sz Giidel 2003). The 
results of these investigations arc in full agreement, converging to a ~ 2.0 — 2.5 
for M dwarfs. If the power-law flare energy distribution extends by about 1-2 
orders of magnitude below the actual detection limit in the light curves, then 
the entire emission could be explained by stochastic flares. The coronal heating 
process would thus be one solely due to time- dependent heating by flares. 
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Table 1. Stellar radiative flare energy distributions 



Star sample 


Photon energies 


a 


References 


M dwarfs 


0.05-2 keV 


L52±0.08 


Collura et al. (1988) 


M dwarfs 


0.05-2 keV 


1.7±0.1 


Pallavicini et al. (1990) 


RS CVn binaries 


EUV 


1.6 


Osten & Brown (1999) 


Two G dwarfs 


EUV 


2.0 2.2 


Audard et al. (1999) 


F-M dwarfs 


EUV 


1.8 2.3 


Audard et al. (2000) 


Three M dwarfs 


EUV 


2.2-2.7 


Kashyap et al. (2002) 


AD Leo 


EUV&O.l-lO keV 


2.0-2.5 


Giidel et al. (2003b) 


AD Leo 


EUV 


2.3 ±0.1 


Arzner & Giidel (2003) 



3.4. Observables of Stochastic Flaring 

In order for stochastic flaring to be an acceptable coronal heating mechanism, 
a number of observables should be correctly reproduced. Flares develop char- 
acteristically in emission measure (EM) and temperature T (fast rise to peak, 
slow decay). The individual flare histories cannot be resolved in the data but 
the statistical, time-integrated thermal histories of the ensemble of flares reflect 
in the time-averaged differential EM (DEM) distribution. The flares follow a 
relation in which the peak luminosity (or EM, hence the emitted X-ray energy) 
is a function of the peak T (Feldman, Laming, & Doschek 1995). For a simple 
flare shape, integration of this dependence over time and over the flare energy 
distribution yields an analytic expression for the DEM: Below the turnover, 
DEM = ri^dV/dlogT oc T^^/^^ ^^ere tt and r„ are the (e-folding) decay times 
of the flare T and electron density, respectively {tt/tu = 0.5 — 2, Reale et al. 
1993). Above the turnover, apart from small correction terms, DEM oc T~^^'^~°'^ 
where 6 ~ 5 it 1 is derived from the Feldman et al. relation. Such DEMs are 
reminiscent of those found from X-ray spectroscopy (Telleschi et al., these proc). 

Flares are often defined in light curves as excess emission above an envelope 
that itself is defined as "quiescent emission". This definition is problematic 
for two reasons, i) It is strongly dependent on the achieved signal-to-noise 
ratio. Any light curve shows quiescent characteristics at sufficiently small S/N 
or time resolution, ii) The superposition of stochastically occurring flares equally 
produces a lower-envelope emission level that is always present. 

Fig. 3 shows simulations of superimposed flares with a = 2.2, compared 
with observations. The large flare observed on Proxima Centauri (Giidel et al. 
2003c) was used as a shape template. The first example uses flare energies spread 
over a factor of ten only, the second uses three orders of magnitude. In the latter 
case, the curve has smoothed out to an extent that it is dominated by "quiescent" 
emission. The individual peaks are merely the peaks of the most energetic flares 
in the ensemble. The available stellar area may constrain the energy range of 
stochastic flares; in small, late-type M dwarfs, flaring active regions may cover 
a significant fraction of the surface (Reale et al. 2003), limiting the number of 
simultaneous flares (upper example in Fig. 3, e.g., Proxima Centauri). 

Spectroscopic density measurements of a time-integrated, stochastically 
flaring corona should then yield values equivalent to the values derived from 
time-integrating the spectrum over a large flare. This is indeed the case. While 
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Figure 3. X-ray light curves. Left panels: Synthetic hght curves from 
superimposed flares; largest to smallest amplitude = 10 (upper) and 
1000 (lower plot). A flare on Proxima Centauri (Giidel et al. 2002) has 
been used as a template, shown filled in the lower plot for the maximum 
amplitude contributing to the light curve. Right panels: Observed X- 
ray light curves for comparison. Proxima Centauri (upper, Giidel et 
al. 2003c) and YY Gem (lower plot, including eclipse; after GOla). 

large stellar X-ray flares achieve peak electron densities of several times lO^"*^ cm~^ 
the time-integrated X-ray spectrum of the Proxima Centauri flare described in 
Giidel et al. (2002, 2003c) yields a characteristic density of logng ~ 10.5 ± 0.25 
derived from O VII, which compares favorably with densities found for magnet- 
ically active main-sequence stars (Ness et al. 2001, 2002). 



4. Conclusions 

Magnetically active stars frequently show bi-modal coronal structure: i) Compact 
"active regions" at small heights appear to be "solar-like" but often require pres- 
sures and densities comparable to solar flare pressures/densities; some compact 
sources have explicitly been seen to develop into expanding flares in radio VLBI 
observations, ii) Very extended, globally well ordered magnetic structure seen in 
X-ray eclipses (as hot gas) or in radio VLBI (as a non-thermal electron popula- 
tion) reaches size scales of several times the stellar radius. These structures may 
be related to cooling, expanding "flare remnants" or large, van Allen-Belt like 
closed magnetospheric structures. Their thermal, density, and temporal prop- 
erties are unlike typical structures in the non-flaring Sun but are reminiscent of 
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flaring plasma. We have discussed expected properties of a stochastically flaring 
corona, such as light curves, the DEM, and the characteristic densities. While 
these predicted observables deviate considerably from values of the non-flaring 
Sun, they are reminiscent of observations of magnetically active stars. 

Why should enhanced flaring occur on magnetically active stars, or in other 
words, why should stars with an enhanced level of magnetic field production 
be prone to develop flare- like magnetic coronae? Wc speculate that, as the 
magnetic dynamo increases toward more active stars, the surface filling factor 
of strong magnetic fields increases, increasing also the volume filling factor of 
closed loops in the corona. As the magnetic loops arc packed closer together 
than on the Sun, they interact more frequently, producing a higher rate of field 
annihilation and reconnection, thus generating flares at a higher rate. These 
energy release events liberate magnetic energy and heat chromospheric plasma, 
evaporating it by overpressure into the magnetic loops. One thus expects both a 
larger average density in the corona and a higher temperature if the flare rate is 
higher. Increased densities and temperatures both are, in this picture, a natural 
consequence of a denser packing of magnetic fields in the corona, induced by 
a more efficient dynamo in stars with more rapid rotation. The determining 
factor of the coronal temperature and emission measure (hence the luminosity) 
is thus the flare frequency (above a certain energy threshold), which in turn is 
a consequence of the amount of magnetic flux production in the stellar interior. 

Acknowledgments. I thank Philippe Andre, Arnold Benz, Robert Mutel, 
and Marek Siarkowski for providing figures reproduced in this paper. 
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